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Abstract 

Contribution of simplest loops is studied in the Local Reggeon Field Theory with a super- 
critical pomeron in the nuclear matter. It is shown that the latter transforms the supercritical 
pomeron into the subcritical. Renormalization of the intercept and difficulties related to the fact 
that the renormalized intercept is complex are discussed. Numerical results with the conventional 
parameters are reported. 

1 Introduction 

With the advent of QCD much effort has been applied to study high-energy hadron-nucleus scatter- 
ing within its framework. However QCD can reliably describe only the hard region of the dynamics. 
Soft processes, contributing to the bulk of the total cross-section, are more difficult to treat. The 
most sophisticated approach is to study soft processes within the perturbative QCD approach to 
small x phenomena, based on the Balitsky-Kovchegov equation, which sums all fan diagrams with 
self- interacting BFKL pomerons [H [21 [3] . However this approach is based on several approxima- 
tions: a large number of colours N c , fixed small QCD coupling constant a s and, most seriously, 
neglect of loop diagrams. This latter approximation can be justified if the parameter 7 = a s exp Ay 
where y is the rapidity and A the pomeron intercept is small. Then for a large nuclear target, such 
that .A 1 / 3 7 ~ 1, the tree diagrams indeed give the dominant contribution and loops may be dropped. 
However with the growth of y the loop contribution becomes not small and although tree diagrams 
are still relatively enhanced by factor A 1 / 3 , disregard of loops cannot be rigorously justified. 

Calculation of the loop contribution within the perturbative QCD appears to be a formidable 
task. So it seems to be worthwhile to start with a simpler approach, that of the old local Reggeon 
Field Theory (LRFT) with a supercritical pomeron. Such a study, apart from its possible lessons for 
the modern QCD approach, also has an independent value, since the old LRFT with phenomeno- 
logical parameters describes not so badly the soft dynamics of high-energy strong interactions. In 
fact in many respects it does it better that the perturbative QCD, which encounters severe diffi- 
culties due to a very large value of the BFKL intercept, hardly compatible with the experimental 
data. Also in LRFT applications to hadron-nucleus [4] and nucleus-nucleus [5] collisions were 
restricted to the approximation, in which the contribution from loop diagrams is neglected. If 
additionally the slope of the pomeron a' is taken to be zero, then the theory is effectively living 
in zero transverse dimensions and allows for the explicit analytic solution. Still remaining in zero 
transverse dimension the influence of loop diagrams has also been studied, both theoretically old 
ago in [HI El E] and numerically recently in [9j . This influence turns out to be decisive for the 
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asymptotic behaviour at large energies and transforms the supercritical pomeron into a weakly 
subcritical one with the effective intercept oc — exp(l/A 2 ) where A is the triple-pomeron coupling 
assumed to be small. 

Unfortunately generalization of these beautiful results to the realistic case of two transverse 
dimensions is prohibitively difficult. To start with, one is forced to introduce a non-zero value of 
the slope: otherwise the loop contribution is divergent in the impact parameter. However even in 
the tree approximation solution of the model with a' ^ is only possible numerically. Second, the 
model in cIt = 2 needs renormalization in the ultraviolet. And, most important, the method used 
to solve the model in dj- = 0, which is to study the corresponding quantum-mechanical system and 
the relevant Schroedinger equation, is inapplicable in the realistic case, since instead of ordinary 
differential equations one arrives at equations with functional derivatives. In fact summation of all 
loop contributions is equivalent to a complete solution of the corresponding quantum field theory, 
the task which seems to be beyond our present possibilities. So at most one can hope to obtain 
some partial results which might shed light onto the properties of the model with loops. There were 
several attempts to study the high-energy behaviour of the LRFT with a supercritical pomeron 
using different approximate techniques and giving contradicting results. In [TO] it was claimed that 
a phase transition occurs at all values of the renormalized intercept = a^(0) — 1 > which 
leads to a theory which violates the projectile-target symmetry and so physically unacceptable. 
On the other hand in [11] it was found that the phase transition takes place only at values of 
e ( r ) greater than some critical value . At < ei the theory essentially corresponds to 
that with a subcritical pomeron and cross-sections vanishing at high energies. At > eP cross- 
sections grow as log 2 s. No violation of the projectile-target was found. These results were obtained 
by approximating the transverse plane as a two-dimensional lattice and starting from the known 
asymptotical solution at each lattice site without intersite interaction. However in this approach 
transition to the continuous plane requires knowing the single-site solutions at large values of the 
triple pomeron interaction constant, which are not known. 

In the present study we take a different approach. Instead of trying to solve the model for 
the purely hadronic scattering we consider the hadron-nucleus scattering and propagation of the 
pomeron inside the heavy nucleus target. Moreover to avoid using numerical solution of the tree 
diagrams contribution with diffusion in the impact parameter, we concentrate on the case of a 
constant nuclear density which allows to start with the known analytical solutions. The point which 
is decisive for the following derivation is that the nuclear surrounding transforms the pomeron from 
the supercritical one with intercept e > to a subcritical one with the intercept — e. Then Regge 
cuts, corresponding to loop diagrams, start at branch points located to the left of the pomeron 
pole and their contribution is subdominant at high energies. As a result the theory acquires the 
properties similar to the standard LRFT with a subcritical pomeron and allows for application of 
the perturbation theory. For a finite nucleus it leads to cross-sections with tend to a constant value 
at large energies. For the physical nucleus with a long distance tail it seems to lead to cross-sections 
growing as log 2 s in agreement with [IT] . 

In fact we cannot prove all these statements rigorously, since we are not able to sum all the 
loops but only a certain simplest subset of them. However the results we find seem to be rather 
general. 

The paper is organized as follows. In the next section we introduce the model and remind some 
relevant known results for the zero-slope case a' = 0. In Section 3 we reformulate the model to 
describe the loops in the nuclear surrounding. Then we calculate the contribution from simplest 
loops and discuss their renormalization in Section 4. Next in Sections 5 and 6 we sum insertions 
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of any number of simplest loops into the scattering amplitude and Green function by solving the 
relevant Dyson equations. Section 7 gives some numerical illustration of the loop influence for the 
amplitude and Green function with physically chosen parameters of the theory. Finally Section 8 
draws some conclusions. 



2 The model. Zero-slope results 

We study the LRFT model based on two pomeron fields 4>(y,b) and cjy(y,b) depending on the 
rapidity y and impact parameter b, with a Lagrangian density 

L = 0^0 + A0*0(0 + f ) + g P <p. (1) 

Here in the free part 

S=\0 y -K ! )+a'Vl + e, (2) 

where e is the intercept minus unity and a' is the slope. The source term describing interaction 
with the nuclear target is 

p = AT(b)5(y). (3) 
It is assumed that for a supercritical pomeron e > and A < 0. The action is 

A = J dy J d 2 bL(y,b) (4) 

and the generating functional is 

Z = f D</>D</Je A . (5) 



The free Green function in the momentum space is 

G^\y,k) = 0(y)e^~ a ' k2 \ (6) 
The classical equation of motion are 

— = -d y </J + a'Vltf + e$ + Xcfi 2 + 2A# f + gp = (7) 

00 

and 

j-^ = d y cP + a'V 2 b (t> + e0 + A0 2 + 2A0 f = 0. (8) 
From the latter equation we find = and the equation for 0^ takes the form 

d y 0t = a'V^0t + e0t + A0t 2 , (9) 

with an initial condition 

0t(y = 0) = gAT(b). (10) 

Equation ([9]) describes evolution of the pomeron field in rapidity and its diffusion in the impact 
parameter inside the nucleus. In the approximation a' = equation for 0^ simplifies to 

d y ^ = e0t + A0t 2 , (11) 
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which is trivially solved for each given b 

1 - \gAT{b)\[ety - 1 



#'<»,»)- 9 -^^ r =«(*,»)■ (12) 



The scattering matrix is 

T(y,b) = gi{y,b). (13) 

3 Loops in the nuclear surrounding 

3.1 Transformation to the nuclear background 

In this section we shall analyze the structure of the model beyond the tree approximation in the 
nuclear background. To locate loops we make a shift in field (fy: 

4>\y,b) = 4(y,b)+^y,b) (14) 

and reinterpret our theory in terms of fields <p and <f>\. 
The Lagrangian density becomes 

L = (4 + £)S4> + X(4 + 2 4> + W 2 (4 +0+ 9P<t>- (15) 
Terms linear in <f> vanish due to equation of motion for £. We are left with 

L = <j)\ (S + + A£0 2 + + 4>) ■ (16) 

This Lagrangian corresponds to a theory in the vacuum with the pomeron propagator in the external 
field/(6,y) = 2A£(y,&) 

G ( f 0) = -(S + 2A£)-\ (17) 

the standard triple interaction and an extra interaction described by the term A£0 2 . This new 
interaction corresponds to the transition of a pair of pomerons into the vacuum at point (y, b) with 
a vertex A£(y, b), see Fig. 1. 

Loops in the Green function may be formed both by the standard interaction and the new one. 
In the latter case they are to be accompanied by at least a pair of of standard interactions. Diagrams 
with a few simple loops in the Green function are illustrated in Fig. 2. One immediately observes 
that a loop formed by the standard interaction has the order A 2 /q' and requires renormalization. 
A loop formed by a new interaction has the order A 3 /a' and is finite. 

The amplitude is obtained as a tadpole g < 0|(y,6) >. The simplest diagrams for it contains 
one loop and are shown in Fig. 3 a, b. Diagrams with more loops are shown in Figs. 3 c, d. 

3.2 The Green function in the external field 2A0^(y) 

In the construction of loops there appears a Green function G^\y,b\y'b') in which the pomeron 
fields are coupled to the external field 

f(y,b) = 2Xav,b). (18) 
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This Green function satisfies the equation 

(e + a'Vl)Gf (y, b\y' , b') + f(y, b)Gf (y, b\y' , b') (19) 



dGf\y,b\y'b')_ (0) {0) 



dy 

with the boundary conditions 

Gf (y,b\y',b')=0, y-y'<0, G^(y\b\y\b') = 5\b -b'). (20) 

In the general case the Green function Gf can only be calculated numerically, just as the 
external field f(y, b). Its analytic form can be found in two cases. 
If the slope a' = then obviously we have 

Gf\y,b\y>,b>) = 5 2 (b-b')Gf\y,y',b), (21) 
where Gj°\y,y',b) is the Green function in the zero-dimensional world which satisfies 

= eG f { y, y>, b) + f{ y, b)G f { y, y>, b) ( 22 ) 

with the boundary conditions 

Gf\y,y',b) = 0, y-y'<0, Gf\y',y',b) = l (23) 

(o 
/ 

Gf\y,y',b) = e J>^^ b » = e-<y~y') ( ^Z^fZ Y = e^'^, (24) 



and f{y,b) given by (|18p and (|12p . Solution of (|22p for Gf is trivial. It is easy to find that 

a-(a-l)e-^'\ 2 , (v _ v <)P 2 (y^ 



a — (a — l)e e ^ y P 2 (y) ' 

where 

. = (25) 

e 

and 

p( y ) = l + a (e £!/ -l). (26) 

It is remarkable that at large y the Green function Gj°\y,y',b) ~ exp(— ey), that is behaves as the 
free Green function with the opposite sign of e and so vanishes at y — > oo. In contrast to the free 
Green function, it corresponds to a subcritical pomeron. 

The second case which admits an analytic solution for the Green function is that of the nuclear 
matter, that is the case when the profile function is constant in all transverse space 

T(b) = T . (27) 

Physically this case corresponds to the behaviour of the Green function at values of of the impact 
parameter well inside the nucleus, where the variation of T{b) is small. Then the classical field £ 
and the external field / = 2A£ become 6-independent and the equation for the Green function in 
this field takes the form 



^^'^ = ( e + a'VDGf (y, b\y', b') + f{y) G f (y, b\y', b') (28) 
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with the same initial condition (J20J). This equation can also be easily solved. We present 



G 



(y, 61^, 6') = - z/', & ? 60e-^' rfyi(e+/(s/l):) = - 6, (y, Z/', &), (29) 



where G^ (y, y', b) is the Green function in zero-dimensional world (|24p satisfying G^ (?/, y', b) = 1 
Then we obtain an equation for X: 

d y X(y,b,b') = a'V 2 b X f (y,b,b') (30) 

with an initial condition 

X(p,b,tf) = 5 2 {b-b') (31) 

and thus with the solution 

X(y-y',b,b') = , j± e-^ 2 /(*<*' (v-v'». (32) 
VKCt\y — y') 

So in the nuclear matter the Green function in the external field takes a simple form 



%',&') = i e -<y-v')-Q>-y?/^{y-y')) ( a Z (° Z ^ ] , (33 ) 

where we used (|24p . In the momentum space we find 

Gf\y,y',k) = e -(y-y')(^ 2 ) ( a Z ( a Z 1 Jf^ V = e (»')(^ 2 )^. (34) 
1 \a — (a — l)e~ e y J p (y) 

Note an especially simple case when a = 1 and the Green function in the external field coincides 
with the free Green function with the opposite sign of e: 



Gf{y,y',k) =e-U-«^ + «*">. (35) 

1 a=l 

Then the theory formally corresponds to a subcritical pomeron model with an additional interaction 
shown in Fig. [1] 

4 Lowest order loops 

4.1 The lowest order loop in absence of the nucleus and its renormalization 

Let us study the simplest contribution from loops, the single loop in the Green function. In the 
second order in A we get for the Green function 

G^(y,k) = -2X 2 J d yi dy 2 ^G^H y -y 1 ,k)G^H yi - y 2,h)G^( yi - y 2,k-k 1 )G^( y2 ,k). (36) 

The minus sign comes from the fact that the simple Green function is in fact -G(°\ so that 111 
the second order we get * )(-G<°>) = -G( 2 \ Factor 2 comes from the contraction of 

(jy (yi, b±)(j) 2 (y2, 62)- Integration over k% is done trivially 

E(y,*) = -2A 2 / pl G W(yM)G(°Hy,k- h) = -C-e 2 ^ 2 ^ 2 , (37) 
J {2tt) 2 y 
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where 

c = At- ( 38 ) 

Ana' 

Yi(y,k) is finite as it stands but leads to divergence at y = when substituted into f|36|) . So it 
requires renormalization, which is achieved by introducing an extra mass term Aec/rip into the 
Lagrangian. The renormalized self-mass is therefore 

S«(y,fc) = E(y,fc) + Ae5(y) (39) 

and we have to require that 

f V dy'S (r) (y', k) = finite terms - C In — + Ae < cx), (40) 

J Drain 

So we conclude that 

Ae = -C]n(c R y min ), (41) 

where y m j n is the cutoff at small values of y and c R is an arbitrary finite constant. Thus we find 
that the renormalized self mass is 



tf r \y,k) = -c(i e ^-(V2Xs/fc 2 +5(y)Hc R y mm )). (42) 

\y J 

or, in the impact parameter space, 

rt r \y,b) = -C(^^e 2 ^~ b2 ^ a ^+5(y)5 2 (b)Hc R y mm )). (43) 

The Dyson equation for the Green function becomes 

G(y,0,k)=G(°\y,0,k)+ dy x G^\y - m ,k) dy 2 ^ r \vi - y 2 , k)G(y 2 , k). (44) 

Jo Jo 

It is trivially solved by the Laplace transformation. We introduce Laplace transforms 

G(E,k)= / dye-BvGiyAk) (45) 
Jo 

and similarly G^(E, k) and Y;( r \E, k). In terms of the Laplace transforms 

G(E, k) = Jj (46) 



From ([6]) and (l42j) we have 



^ E ^= E- t + ^ - (47> 



and 



Vmin y 

C{ln(E -2e + a'k 2 /2) + C E -lnc R ). (48) 
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(Note that terms with ln(y m j n ) cancel). So the Laplace transform for the Green function is 

G-\E, k)=E-e + a'k 2 - C[ln(E - 2e + a'k 2 /2) + C E - lnc R ]. (49) 

Its singularities in the complex E'-plane consist of a branchpoint at E = 2e — a'k 2 /2 and a pole 
Eo(k) determined by the equation: 

G- 1 (E ,k) = (50) 
At small A Eq. ([50]) gives two complex conjugate poles 

E^ = e - a'k 2 + C(ln(|e + a'k 2 /2\) + C E -lnc R ± in). (51) 

These poles appear on the physical sheet due to the abnormal sign of the self-mass contribution, 
which is a consequence of the imaginary coupling in the original Lagrangian of the LRFT. The 
asymptotic of the Green function is ~ exp(2/xy) in accordance of the contribution of the Regge cut 
for a supercritical pomeron. According to (|4ip the unrenormalized pomeron intercept (depending 
on the cutoff) is 

eo(yo) = e - C\n{c R y min ). (52) 

It goes to infinity as y m in - » 0, since C > 0. As to the renormalized intercept its definition may be 
chosen in different ways. Traditionally it may be taken as the position of the pole in E in the Green 
function at k = 0. Then it is given by (|50|) . complex and c^-dependent. One observes the difficulty 
in doing renormalization in the standard way by requiring the zero order propagator to have the 
pole at a chosen energy Eo(k = 0): unless Eq > 2e there are two complex conjugate values of Eq 
and if one chooses Eq > 2e contributions from higher order self mass will again make it complex. 
So the structure of perturbative singularities prohibits choosing the zero-order propagator to carry 
the pole singularity, since there are two. However from the physical point of view such a choice is 
not necessary Inclusion of more pomeron exchanges shifts the dominating singularity further to the 
right, so that the simple pole contribution to any physical process becomes completely obliterated. 
It is more reasonable to relate the value of the renormalized intercept (and thus fix the arbitrary 
constant c R ) to some physical observables. In the following we shall see that in the scattering on 
the nucleus there appears a possibility to define the renormalized intercept from experimental data 
at high energies, with a certain choice of c R . 

4.2 Lowest order loop in the nuclear matter 

The contribution from the lowest order loop (Fig. 2 a) in the nuclear matter to the Green function 
in is 

Gf(y,0,k) = J dy 1 dy 2 Gf\y,y 1 ,k)J: 1 (y 1 , y2 ,k)Gf ) (y 2 ,0,k). (53) 
where the Green functions G^p are given by (1241) and 

Xi( yi ,y 2 ,k) = -C(^—e- 2 ^-^-( 1 / 2 ^-^ a ' k \ 2j y^ (54) 

K yi-y2 ' 

is the renormalized 2nd order self mass in the nuclear background. The subtraction term is easily 
calculated: 

Gfsubtr = -CHcRVmm) J dy 1 Gf(y,y 1 ,k)Gf ) (y 1 ,0,k) = -CHc R y min )yGf (y,0,k). (55) 
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The main term is given by 

-CGf(y,0,k) /V [ Vl dy 2 ^—e^^ 2 ^. (56) 
1 Jo Jo 2/1-2/2 P{yi) 

Changing to integration variable z = y\ — y 2 we find 

Gfl ain {y,0,k) = -CGf\y,0,k) J o y ^W(z), (57) 



where 



W (z) = e^\^ 2 ) I" dyi P {V 9 ) ~ Z) ^ e^' k2 W(z). (58) 
Jz p 2 {yi) 



We note that W(0) = y so that we can write 

Gf mam {y, 0, k) = -CGf\y, 0, k){ I" -(W(z) - y )+yln^-}. (59) 

The second term just substitutes in the subtraction contribution y m in by y. The first term can be 
integrated numerically since it converges at z = 0. 

Function W(z) can be calculated analytically. Passing to integration variable u = e eyi and 
doing simple integrals we find 

vtrt \ -<*( i gf ~ I 2 i n 2("i+g) (gl ~ g) 2 «2 - m \ 

W(z) = e [y — z-\ 5 — In — ; r -. I. (60) 

V e q* u\{U2 + q) eq ^ 2 _|_ q)(u\ + q) > 

where u\ = e ez , u 2 = e €y , q = 1/a — 1 and q\ = qu\. 

Note that after renormalization the 2nd order contribution to the Green function is also finite 
in the nuclear background. 



4.3 3d order loop Fig. 2 b 

The new loop shown in Fig. 2 b, which appears in the nucleus, has formally order A 3 . However it 
is proportional to AT (Jo), so that the extra A combines into the characteristic parameter o, which 
is not small. So in fact it has the same order of magnitude as the simple loop of Fig. 2 a. 
The contribution to the pomeron self-mass corresponding to Fig. 2 b is 

£ 2 (j/2,y3,A0 =^ d Vl J ^Gf\y 2 ,y 3 ,k-k l )Gf\y 2 ,y 1 ,k 1 )Gf\y s ,y 1 ,k l )^y 1 ). (61) 

Here it is assumed that y 2 > j/3. The final factor 4 combines factor 4 from contractions of 4> 2 {yi) 
with <jy (y 2 ) and <fi^ (2/3), factor 2 coming from interchange 2^3 and factor 1/2 left from 1/3! 
in the development of action in the third order in interaction. The sign " +" comes from the odd 
number of propagators in the contribution to the full propagator with a minus sign. Using explicit 
expressions for the propagators and £ and definition (|25p we rewrite (|6ip as 

Z 2 (y 2 ,y 3 ,k) = -4\ 2 ae f" d^^M^-i/O-te-MWfc 2 / ^ e -2a'fc?fc-y 1 ) + 2(, 2 -, 3 K^i. 

A) V KV2) J (2vr)^ 

(62) 
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The integral over k\ is 

1 / fc 2 (V2-V3) 2 

g 2 (V2-Vl) 

8Tra'(y 2 -yi) 

so that in the end we obtain 

£^,,3,*, = -pL 1 «^-*.^ f _^_/ (si)e -<». + «'' 2 fe5 (63 ) 
27ra'p 4 (y 2 ) Jo 2/2 - 2/1 

It is trivial to see that £2 is finite and does not need any more renormalization. 

5 Amplitude in the nuclear matter approximation 

5.1 Lowest order 

We start from the lowest order approximation, when at fixed b the scattering amplitude is 

T^(y,b)=g 2 AT(b)^- (64) 

The total forward scattering amplitude is obtained after integration over all b: 

T(°)(y)= / d 2 bT^(y,b). (65) 



T {0) (y,b)= I A* e «f(o) (yifc)> (66) 



This expression can be rewritten in the momentum space, once we introduce the Fourier transform 
of T(y, b) by 

d 2 k 
(2nf 
as 

T^(V) = J ^e ikb fW(y,k) = r<%,0). (67) 
In the nuclear matter T(y, b) does not depend on b: 

T^(y,b)=g 2 AT - pe ' ! 



p(y) 

so that its Fourier transform is 



f^(y,k) = {2^) 2 5 2 {k)g 2 AT Q - 



p(y) 

From (|65p the forward scattering amplitude is 

r(%) = / d 2 bt(y) = itRWAT^ = g 2 A^- y (70) 

where we used that To = \/{itR 2 A ). Due to (fBTl) this of course implies that {2 r K) 2 5 2 {k = 0) — > ttR 2 a . 
In the limit y — ► 00 the lowest order amplitude tends to a finite value 

^(|/U = ^ = ^|. (71) 
a |A| 
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5.2 2nd order 

In the next order we have 



fW(y,k)= / dy 1 dy 2 Gf ) (y,y 1 ,k)E 1 (y 1 ,y 2 ,k)f^(y 2 ,k) 



= (2TT) 2 5 2 (k)g 2 AT J dy 1 dy 2 Gf\y,y 1 ,0)E 1 (y 1 ,y 2 ,0) p 
so that the forward scattering amplitude is 



/ 



T^(y) = TrR 2 A J dy 1 dy 2 G f (y,y 1 ,0)^ 1 (y 1 ,y 2 ,0)T^(y 2 ,b). 
The subtraction term in T,^ gives a contribution 

T sll(y) = -CirR 2 A Mc R ymin) £ dy 1 Gf ) (y,y 1 ,0)T^(y 1 ,b). 
Using ((SH) and ((24]) we get 

T sub(v) = -Cg 2 AGf{yM\n(c R y mm )J\y 1 p(y l ). 



We define 



Then we find 



I(y, z) = J" dy lP ( yi ) = a -(e^ - e* z ) + {y-z)(l-a). 

T sub(y) = -Cg 2 A^ln(c R y min )I(y,0). 
The main term is 



(2) 2 e f V f Vl 

T rnain(y) = ~ C 9 A -J7-T / / dy 2 e e 

P (y) Jo Jo 



yi -V2P 2 (yi)' 

Passing to integration variables yi and z = y\ — y 2 we find 
where 

fj/ p 3 (?/i-^) 



*(z)=e« fW 

J 2 



p 2 (yi) 



Note that 

#(0)=I(y,0). 

Thus we find 

e ey r /"^ dz 



T^(y) = -^ 2 ^{ /; t(^) - *(»)) + m jiL/fc, 0)} 

The second term just changes y m i n to y in the subtraction term. 
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Function ^f(z) can be calculated analytically in the same way as W(z) (see Eq. (|58p ). One 
finds 

lTf / \ KVi*) i o ux-l (ui-1) 2 U2 {u 1 +q) 

+ g q ^ Ul - 1) \ Ul ' U2 r . (83) 



where as before u\ = e ez , u 2 = e ey , q = 1/a — 1 and q\ = qui. The remaining integration over z 
has to be done numerically. 

It is not difficult to find the asymptotic behaviour of T^ 2 \y) at large y. At y — > oo, asymptoti- 
cally 

p(y) = ae*y, 7(y,0) = -e e « (84) 

e 

so that in this limit 

a 2 A 

T {2) (y)sub = -C y — Hc R y min ) (85) 
ae 

The asymptotic of the main term comes from large values of z inside the integral. At large z 

$>( z ) _^(o) = -e~ 2ez ) (86) 

and so 

|"^(*(z)-*(0)) ~ V«(Ei(-2e7/)-ln(2ej/)-C £ ) ~ - V"(ln(2ey) + C B ). (87) 
Thus asymptotically 



?i 2 L(y) = C^-(ln(2 e2/mm ) + C £ ). 



ae 

In the sum the terms with lny. m j n cancel and we find that in the limit y — > oo the loop 
contribution also tends to a constant 

T {i \y) y ^ 00 = -C 9 ^(\n C f-C E ). (89) 
ae V 2e / 

The ratio first to second order is at y — > oo 

'"'W-^IU-^Cnf-C,). (00) 
Remarkably it does not depend on A. 
5.3 3d order 

In the third order the forward scattering amplitude is 

T^(y) = ttR 2 a J dy 2 dy 3 G f (y,y 2l 0)Z 2 (y 2 ,y 3 ,0)T^(y 3 ,b). (91) 

Using (|63p at k = and interchanging the order of integration in y 2 and internal integration in y\ 
inside £2 we find 

T^(y) = -g 2 A^4- dy,^- -J^ p ^y l)e -^ dy * (92) 
2ira'p z (y)Jo P z {y2) Jo 2/2 - yi Vyi PUW 
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It is not difficult to find the asymptotic of T^ 3 \y) at large y. We present (|92p in the form 

T^(y) = -g^A^^T- I" dV2-ff^F(y 2 ), (93) 
2ira f p 2 (y) Jo p 2 {y 2 ) 



where taking y 1 = y 2 (3i and y 3 = y 2 (3- 



3 



F(y 2 ) = V2 C T^W(y 2 /?i)e-^ A /' 
J o 1 — Pi Jpi 



We are interested in the behaviour of F(y 2 ) as y 2 >> 1, since this governs the behaviour of the 
integral over y 2 in ([93j) . As y 2 — » oo we have p(y 2 (3) — ► aexp(ey2/?) so that 

= a a W jT 1 ^L-e 2 ^ j£ d/3 3 = ^ (e 2 ^ - l) . (95) 



Putting this asymptotic in (|93p we obtain 

T»( ! /)»» 1 = -fA^e-' £ <ly 2 e™ = (96) 
Recalling that at y >> 1 T(°) = ff 2 A/a we find the ratio 

-~ m 



r«%) 



y»l A'Ka'e 



Note that in the particular case a = 1, T^ 3 ) can easily be found explicitly at all y. In fact trivial 
integrations give 

TS(y) = V^(l-e--) 2 . (98) 
The full ratio T/T^ from both loops, Fig. 2 a and b, turns out to be 

r(y) = (y) + (y) = ( In ^ + 1 - C E ) . (99) 

47ra e v 2e / 

The bracket is universal. So one may define the renormalized self-mass by requiring that this 
bracket is zero, which implies that the loop correction vanishes at large y: 

In — - C E + 1 = 0, or c R = 2ee c ' B " 1 . (100) 

This allows to experimentally determine a and hence e from the behaviour of the amplitude at large 
y. In this way we may define the renormalized intercept in a physically reasonable manner. 

Note that with this choice of cr the renormalized intercept formally defined as the position 
of the pole of the propagator in the vacuum at k = is given by the solution of the equation (|50p 
which reads 

i n - 2e 

e (r) - e - C[ln + 1] = (101) 

At small A it is complex and close to e. 
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5.4 Random phase approximation 

One may try to sum all primitive loop insertions into the amplitude ('random phase approxima- 
tion'). In this approximation the amplitude T(y) in the nuclear matter satisfies the Dyson equation 

T(y) = r(°)(y) + £ dy x jT dy 2 Gf{y, y x , 0) (x x (y x ,y 2 ,0) + S 2 (l/i,!fe, oj)T(y 2 ), (102) 

where £i and E 2 are the 2nd and 3d order self-mass contributions studied above. 
We present 

T{y)=T^{y)r{y). (103) 

The equation for r(y) reads 

r(y) = l + X x (y) + X 2 (y), (104) 

where X x and X 2 are parts coming from Si and £ 2 respectively. 
We find 

c [y c fy ry 1 i i )P 3 (y2) 

x l(y) = r^ln(c fl y miri ) / dy x p{y 1 )r(y 1 ) — / dy x / dy 2 e e(Vl y2> M y 2 ). 

p{y) Jo p(y) Jo Jo 2/1-2/2 p 2 m) 

(105) 

Again we pass to variable z = y x — y 2 to rewrite the last term as 

C l y -^z) = -^[ y ^Uy,z)-u J { y ^))-^My^)f y * (106) 
Jo z p(y) Jo z V > p{y) J Vmin z 



p{y) 

where 



u(y,z) = e a I* dy /^f. *K (y x -z). (107) 

Jz p 2 {yi) 

Obviously 

wfc/,0) = f V dy x p(y x )r(y x ), (108) 
Jo 

so that the last term in (|106|) changes 2/min. in (|105|) to y. As a result we find that the part X x (y) 
in the equation for r takes the form 

c fy c fy dz / \ 

Xl = ~p(v) Jo d2/ip(yi)r(yi) ~pJy)J Q ~ z) " 0) ) • (109) 
The part X 2 (y) is 

2aeC fy e eyi fV 1 dv* □ e ey2 

X 2 (y) = --^p- / dj/i-^ / ^^P 3 (y 3 )e-^ / ^2/ 2 ^r(y 2 ). (110) 

p{y) Jo p\y\) Jo 2/1-2/3 J yz v\v*) 

In the case a = 1 and neglecting X 2 we can solve the Dyson equation passing to the Laplace 
transforms in rapidity. We find 

T (E) = / _^- 1 , V LL —T. , (HI) 



where (E, 0) and Hp (E, 0) have the same form as (|47j) and (|48p with the opposite sign of e 
and 

T^{E)=g 2 A^. (112) 



J5 
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So the Laplace transform for the amplitude is 

t(f\ = 9 — E±l ni^ 

1 ' E E + e-C[\n{E + 2e) + C E -\nc R \ K ' 

From the expression for T{E) we immediately conclude that in the limit y — > oo the amplitude 
tends to 

T(y)\ y ^ 00 = g 2 AZ, (114) 

where 

Z= e-C[ln(2e) + C E -lnc R y (U5) 

Factor Z — 1 gives the correction from loops. We can define the renormalization condition by 
requiring that Z = 1, which implies for this case (X 2 = 0) 

ln(2e) + C E -lncfl = 0, (116) 

With this choice 

^(E,k) = Cln E + 2e l a ' k2/2 (117) 

and 



T(E) = 9 -^ E + e E+2 . (118) 

6 Pomeron Green function in the random phase approximation 

In this approximation the full Green function in the nuclear matter Gf(y, 0, k) satisfies the Dyson 
equation 

G f (y, 0, k) = Gf\y, 0, k) + j" d Vl dy 2 Gf{y, y ll k)(z 1 (y 1 ,y 2 , k) + E 2 (y 1>y2 , k))G f {y 2 , 0, k) 

(119) 

We present 

G f (y,0,k) = Gf\y,0,k)R(y,k) (120) 

to obtain an equation for R 

R(y,k) = l + Y 1 +Y 2 . (121) 

Here Y\ comes from Ex: 

Y 1 ( y ,k) = -Cln(c R y min ) f d yi R(y lt k) - C f d Vl dy 2 ^—e^'y^^}-R(y 2 ,k) (122) 

Jo Jo Jo 2/1-2/2 V (2/1) 

and j3 = e + a'k 2 /2. Passing to integration variable z = y\ — y 2 we rewrite the last term in Eq. 
CE22D as 

ry A z ry A 7 , . ry A 7 

-c -n( y ,z) = -c -(n(y,z)-n( v ,oj)-cn(y,o) - (123) 

Jo z Jo z V 1 J Vmin z 

where 

n(y,z) = e^ f V dy / { f~. z) R( yi -z,k). (124) 
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Obviously 

Q(y,0) = F d yi R(yi,k), (125) 
J o 

so that the last term in (|123p changes y m in m (|122|) to y. As a result we find 

ry ry A~ / \ 

Yi(y,k) = -Cln(c R y) J d Vl R( yi ,k) - C J — (n(y,z) - J2(y,0)). (126) 
The part Y2 comes from £2: 

Y 2 (y,k) = 

9 ^ Fa &eVl F 1 dy2 3/ n T 1 , e g w /l , 2 (2/1-2/3)2 x 
-2aeC dy! / P (2/2]e y2 / 2 , ^ R{V3,k) exp I -a fc . (127) 

p 2 (yi) Jo 2/1-2/2 -A/2 p 2 (y3) v 2 2/1-2/2 y 

In the general case Eq. (|126p can only be solved numerically. In a particular case a = 1 both 
the Green function and £1 depend only on the rapidity difference: 

p(y) = e ey (128) 

and 

Gf\y h y 2 ,k) = e -(w-w)(^* a ) = ( yi - 2/2, k), (129) 

Si (2/1, 2/2, fc) = -C(—!—e~^-^ +a ' k2 W + _ y2 ) \n(c R y min ) = Si(yi - 2/2, fc). (130) 

They have the same form as in the vacuum with the opposite sign of e. So if one additionally 
neglects the formally 3d-order contribution £2 the Dyson equation f)119f) can be again analytically 
solved by the Laplace transform in the same way as in the vacuum case. Moreover, the solution is 
exactly the same as in the vacuum case with the change e — > — e: 

G~ f l {E,k) = E + e + a'k 2 - C[\n(E + 2e + a'k 2 /2) + C E - lnc R }. (131) 

Singularities of Gf(E, k) in the complex E— plane consist of a left cut starting at E c = —2e — a'k 2 /2 
and corresponding to the standard cut generated by a subcritical pomeron with intercept 1 — e and 
a possible pole whose position depends on the renormalization constant c R and which determines 
the asymptotic at large rapidities. With the choice (|116p corresponding to £2 = we have 

G- f \E, k) = E + e + a'k 2 - Cln E + 2e + a ' k2 / 2 (132) 
J 2e 

The Green function vanishes at y — > 00 as e~ ye ° where for a'k 2 < 2e 

e + a'k 2 < e < 2e + a'k 2 /2. (133) 

If a'k 2 > 2e, the two-pomeron cut moves to the right of the pole, which splits into two complex 
conjugate poles on the physical sheet, and the asymptotic is determined by the cut so that y = 
2e + a'k 2 /2 imaginary parts of the two terms in the denominator are opposite. 

Unfortunately inclusion of £2 spoils invariance with respect to translations in rapidity, and 
solution of the Dyson equation by means of the Laplace transform becomes impossible. 
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7 Numerical illustration 



In this section we report on numerical results for the amplitude T(y) and Green function Gf(y, k) 
which follow from the Dyson equations (|102p and (|119|) with the pomeron self-mass £ = Si + S2 
given by loops of Fig. 2 a and b and realistic values of e, a', A and g For the latter we take the 
standard values, which correspond to the experimental data at comparatively low energies 

e = 0.08, a =0.2 GeV~ 2 , A = -0.48GeV -1 , g = 5.94 GeV^ 1 (134) 

For the transverse nuclear density we choose that in the center of the nucleus with the constant 
nuclear density within a sphere of radius Ra = A 1 /^ 1.15 fm. 

With e, a' and A given by (|134p we obtain from Eq. (jlOip the renormalized values for the 
intercept and slope in the vacuum: 

e W =0.155 ±i 0.139, a {r) = 0.172 t« 0.043 (135) 

Both are complex. The real part of the pomeron trajectory is found to be linear in k 2 with a good 
precision. The imaginary part is not, the value of Im a!^ r ' diminishing from 0.043 at k = to 0.020 
at k = 2 GeV/c. So we are dealing with a supercritical pomeron with a complex intercept and 
slope. 

Next we pass to the amplitude and pomeron Green function in the nuclear background. 

Fig. H] shows the ratio r(y) of the calculated amplitude to the lowest order one as a function 
of y for two values of the atomic number A = 64 and 207. As one observes as rapidity grows the 
ratio initially goes up reaching at its maximum values of the order 2.2 -7- 2.3 depending on A and 
then gradually goes to unity, as determined by the chosen value of cr. The ^-dependence is in 
fact weak, so that the two curves are close to one another. The rather large values of the ratio 
at intermediate rapidities show that for the chosen parameters the loop contribution is not at all 
small. In other words the strength of the triple pomeron interaction is relatively large although it 
effectively goes down at large rapidities. 

In Figs. 07] we show the R(y,k) of the Green function to its lowest order value as a function 
of y for values of k = 0, 1 and 2 GeV/c and A = 64 and 207. For convenience for k = 2 GeV/c 
we plot R\(y, k) = —R(y, k) exp (y(e — a'k 2 /2)j , since for such large k the asymptotic is negative 
and governed by the cut contribution. Again we observe that the ^4-dependence is weak. The 
y-dependence strongly depends on the chosen value of k, which is to be expected, since the cut 
contribution vanishes at k — > 00 much slower, with a twice smaller slope. For all momenta the 
loop contribution is relatively large at all rapidities, which again indicates that with our choice A 
is not small. Oscillations observed in Figs. [5] and [6] illustrate that the 'physical' intercept e( r ) for 
moderate k in fact splits into two complex conjugate values. At larger k the asymptotic is taken 
over by the cut contribution and oscillations disappear. 

8 Conclusions 

We studied the contributions of the two simplest loops in the local reggeon field theory with a 
supercritical pomeron in the nucleus. Our results show that the nuclear surrounding effectively 
transforms the supercritical pomeron with the intercept a(0) — 1 = e > into the subcritical one 
with the intercept — e. As a result, at high energies the pomeron Green function vanishes and 
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contributions from multipomeron exchanges vanish still faster according to the standard predic- 
tions for the subcritical pomeron. With that, for a finite nucleus with a constant profile function, 
the scattering amplitude tends to a constant value at high energies. Contribution from self-mass 
insertions do not change this behaviour, since in the rapidity space they are dominated by the con- 
figuration in which these insertions enter at small rapidities and the bulk of the rapidity is covered 
by the lowest order amplitude. With an appropriate choice of the renormalization constant one 
can make the contribution from loops vanish at high energies. The renormalized pomeron inter- 
cept in this theory (in vacuum) is complex with a positive real part. This makes the conventional 
renormalization technique inapplicable. 

These results have been obtained for a simple case of a constant nuclear profile function T(b) and 
with only lowest order loops taken into account. The limitation to a constant T(b) is technical and 
we believe that the same conclusions can be obtained for a realistic T(b). However this requires using 
non-trivial pomeron Green functions in the external 6-dependent field, which can only be calculated 
numerically. Construction of loops in this case presents a formidable calculational problem, which 
seems to be outside our technical possibilities. The only difference we expect with a realistic Tip) 
is due to its long distance tail. Then at large b > b${y) the amplitude and the Green function 
will no more be damped by the nucleus and become purely perturbative, that is growing as e ey 
at large y. Very crudely we expect bo to be determined by the relation T(bo)e ey ~ 1 which gives 
bo(y) ~ ay where o is an ^4-independent parameter of the dimension of length (a = 0.545 fm with 
the Woods-Saxon nucleus density). Then, again very crudely, the cross-section obtained by the 
integration over all b will become proportional to y 2 in accordance with the conclusions in 
Inclusion of more complicated loops does not present any difficulty of principle either. For small 
enough A their contribution will be small at all rapidities and we do not expect any qualitative 
change. 

Note that with the conventional values (I134p for e, a' and A the loop contribution proves to be 
not at all small, in spite of the rather small value of the effective loop parameter C ~ 0.09. The 
renormalized intercept (|135H obtained as the solution of Eq, (|50p is quite different from its zero- 
order value e = 0.08, which testifies that the coupling constant A is in fact quite large. Because 
of that we cannot claim that with the conventional parameters (1134p our results are complete . 
More complicated loops have to be included to make the results reliable and fit to be used for the 
description of physical observables. 

This is unfortunate, but we consider our results to be interesting mostly from the purely theo- 
retical point of view. They may serve as a basis for treating loop contributions in the perturbative 
QCD. It seems to be advantageous to study them in the nuclear surrounding, which makes the 
high-energy behaviour much more tractable. 

Appendix A: Possible poles of propagators 

Let the denominator in (|49p be 



where we have taken into account that C < 0. We denote 5 = \C\(Ce — utcr) and rewrite (lAlj) as 



Function D(E,k) has a left cut starting at E = 2e — a'k 2 /2, so that it can vanish only to the 
right of the branchpoint. At E — > 2e — a'k 2 /2 from the right D(E,k) — > +oo. At E — > +co also 



D(E,k) = E -e + a'k 2 - \C\ [ln(£ - 2e + a'k 2 /2) + C E - hxc R \, 



(Al) 



D(E, k) = E-e + a'k 2 - 6 - \C\ ln(E - 2e + ak 2 /2). 



(A2) 
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D{E, k) — > +00. The single minimum of D(E, k) as a function of E occurs at 

dD ( E > k ) 1 \ C \ 1 Q (A3) 

that is at 

£™n = e |c| + 2e - ^a'fc 2 . (A4) 

The minimal value of D(E,k) is 

D min (E, k) = e' c ' + e + ^a'A; 2 - 5 - \C\ 2 . (A5) 

Obviously if -D m j n > then D(E, k) does not vanish on the physical sheet of E. If D m i n < there 
are two zeros to the right of the cut. The rightmost zero Eq determines the asymptotic of the 
propagator at large y. The value of D m i n depends on the choice of 5. So if 

5> +e + ^a'k 2 -\C\ 2 (A6) 

the propagator develops two poles to the right of the branchpoint at E = 2e — a'k 2 /2, the larger of 
which is located to the right of E m i n and thus to the right of the unperturbed pole at E = fi — a k 2 . 
This larger pole takes the function of the renormalized intercept. If, on the other hand, relation 
(|A6p is not satisfied, the propagator has no poles on the physical sheet and its asymptotic at y — > 00 
is determined by the branchpoint at E = 2e — a'k 2 /2 

Passing to the nuclear surrounding with a = 1, at k = we find a denominator 

D{E) =E + e- |C|ln^i^ (A7) 
We put E + 2e = 2e£, |C|/(2e) = 77 << 1, so that zeros are determined by the equation for £ 

/(e) = e-r ? ine-i = o. (as) 

The minimum of /(£) occurs at Z, m in = and 

fmin = V ?/ In 77 ~ (A9) 

At 77 < 1 /mm < for 7/ < 0.186683, for which values there exist two zeros 

0<£ 1 <r/<e 2 <i (A10) 

(the limiting value 1/2 is achieved by the larger zero at rj — > 0). This means that poles in the 
propagator exist for small enough r/ and are located between the unperturbed pole at E = — e and 
branchpoint at E = 2e. One can easily see that this result is also true for k 2 > 0. So the poles do 
not change the asymptotic of the amplitude, but diminish the intercept of the propagator, making 
it still more subcritical. 
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Figure captions 

Figure 1: The new vertex for two-pomeron annihilation, which appears after the shift in field (fy . 
Figure 2: Some simple loop diagrams for the pomeron Green function. 

Figure 3: Diagrams with one loop (a, b and two loops (c and d) for the scattering amplitude. 
Figure 4: The ratio r(y) of forward scattering amplitude to its lowest order for Cu and Pb targets. 
Figure 5: The ratio R{y) of the pomeron Green function to its lowest order at k = for Cu and 
Pb targets. 

Figure 6: Same as Fig. [5] at k = 1 GeV/c. 

Figure 7: The scaled ratio R\(y, k) = —R(y, k) exp (y(e — a'k 2 /2)j at k = 2 GeV/c for Cu and Pb 
targets. 
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Figure 1: The new vertex for two-pomeron annihilation, which appears after the shift in field $ 
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Figure 5: The ratio R(y) of the pomeron Green function to its lowest order at k = for Cu and 
Pb targets 
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